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ABSTRACT 

We present spatially-resolved analysis of the temperature and gas density profiles in 6 relaxed galaxy clusters 
at z = 0.4 - 0.54 using long-exposure Chandra observations. We derived the total cluster masses within 
the radius rsoo assuming hydrostatic equilibrium but without assuming isothermality of the intracluster gas. 
Together with a similar study based on the XMM-Newton observations (Kotov & Vikhlinin), we obtained the 
mass and temperature measurements for 13 galaxy clusters at 0.4 < z < 0.7 spanning a temperature interval 
of 3 keV < T < 14keV. The observed evolution of the M - T relation, relative to the low-redshift references 
from the Chandra sample of Vikhlinin et al., follows MscalT^^^ E{zy, where we measure a = 1.02 + 0.20 
and ff = 1.33 + 0.20 for the spectroscopic and gas mass-weighted temperatures, respectively. Both values 
are in agreement with the expected self-similar evolution, a = I. Assuming that the cluster mass for given 
temperature indeed evolves self-similarly, the derived slopes, y, of the high-redshift M - T relation, £'(z)M50() oc 
r^, are y = 1.55 + 0.14 for Tsp^c and y - 1.65 + 0.15 for T,„g. Our results show that both the shape and evolution 
of the cluster M - T relation at z ^ 0.5 is close to predictions of the self-similar theory. 

Subject headings: galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

Scaling relations between the cluster parameters such as 
total mass, average gas temperature, and X-ray luminosity 
can be used for studying galaxy clusters and their cosmolog- 
ical applications (see Voit 2005 for a recent review). Several 
recent Chandra and XMM-Newton studies provide accurate 
measurements for the M - T relation in low -redshift clusters 
JVikhlinin et alJl2005hH Arnaud et al."2005Y While studying 
the scaling relations for low-redshift clusters allows to ad- 
dress many interesting questions, the knowledge of the evo- 
lution of the scaling relations at z > is required in many 
cosmological applications. The X-ray cluster mass measure- 
ments based on application of the hydrostatic equilibrium 
equation (Mathews 1978; Sarazin 1988) rely on observations 
of the cluster temperature and brightness profiles at large 
radii. These are technically challenging observations, espe- 
cially for distant clusters. In the past, the mass measurements 
for distant clusters were often derived assuming T{r) - const 
(e.g., Ettori et al. 2004; Maughan et al. 2005). However, 
this assumption can significantly bias the mass measurements 
JMarkevitch & VikhUnin. 1997: Maughan et aL.2005l 

The first systematic study of the high-redshift cluster M-T 
relation based on the spatially-resolved temperature profiles 
to r * rsoo ('"500 is defined as a radius corresponding mean 
overdensity A - 500 relative to the critical density at a cluster 
re dshift) fro m the XMM-Newton observations was presented 
in lKotov'& Vikhlinin (200j). In this Paper, we expand this 
work by analysing Chandra observations for 6 clusters at 
z - 0.4 - 0.54. All our clusters belong to a sample of lumi- 
nous, dynamically relaxed systems used bv lAUen et alJJ2004h 
for fgas measurements. Statistical accuracy of the observa- 
tions from our sample allows us to trace the surface brightness 
and temperature profiles nearly to rsoo. Chandra's superb an- 
gular resolution minimizes the systematic uncertainties in the 
spatially-resolved analysis. We use the derived gas density 
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and temperature profile measurements to infer the total clus- 
ter masses and the average temperatures. A combination of 
the Chandra and XMM-Newton distant cluster samples pro- 
vides an accurate determination of the shape and evolution of 
the cluster M-T relation at z ^ 0.5. 

All distance-dependent quantities are derived assuming the 
Qm - 0.3, Qa = 0.7 cosmology with the Hubble constant 
Hq - 72km s"' Mpc"'. Statistical uncertainties are quoted at 
68% CL. 

2. OBSERVATIONS AND DATA REDUCTION 

In present analysis, we used the archived Chandra obser- 
vations of z > 0.4 clusters. We selected only those objected 
with the Chandra exposures sufficient for temperature mea- 
surements to r ~ rjoo. We also selected clusters which appear 
dynamically relaxed in their Chandra images. The selected 
clusters are listed in Tabl^ 

Our Chandra data reduction follows the procedure de- 
scribed in Vikhlinin et al. (2005a b). This involves generation 
of the spatially dependent effective area and detector response 
files, modeling of the particle-induced detector background 
following Markevitch et al. (2003), and subtrac tion of the 
residual Galactic foreground as described in .Vikhlinin et alJ 
(2005a). 

All detectable point and small-scale extended sources were 
masked out from the analysis of azimuthally averaged sur- 
face brightness and temperature profiles. The sources were 
detected using a wavelet decomposition technique described 
in Vikhlinin et al] (ll998h . The sources were automatically de- 
tected in the 0.7-2 keV and 2-7 keV energy bands and the 
merged region mask was checked manually. 

All visible large-scale extended X-ray sources were also ex- 
cluded from the data. Such structures are present in 3 clusters. 
In MACSJ0159. 8-0849 there are detectable filamentary struc- 
tures with a size of ~ 2', located at ^ 5' from the cluster center 
(at RA = 30.031°; Dec = -8.854° and RA = 29.965°; Dec = 
-8.921°). MACSJ1621.3+3810 is projected on the outskirts 
of a low-redshift galaxy group (~ 25' from the group cen- 
ter). The X-ray emission associated with the group is clearly 
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TABLE 1 

Best Fit Parameters for Gas Density Profiles and Results of Spectral and Mass Determination 



Cluster 


z 


f!02/«o'' 






a'' 








Pi' 


''500 


^spec 












(kpc) 


(kpc) 








(kpc) 




(kpc) 


(keV) 


(keV) 


(10'" Mo) 


MACS 1423.8+2404. . . 


0.539 


14.9 


28 


1220 


1.86 


0.53 


4.94 


0.02 


1.00 


973 ± 72 


7.02 ± 0.28 


6.27 ± 0.40 


4.56 ± 1.04 


3C295 


0.460 


173.5 


166 


1548 


2.47 


0.63 


4.99 


5.80 


1.00 


840 ± 66 


5.13 ±0.24 


4.52 ±0.34 


2.63 ± 0.66 


MACSJ0159.8-0849 . . 


0.405 


27.3 


247 


2308 


1.84 


0.69 


5.00 


23.74 


1.00 


1324 ± 127 


9.59 ±0.50 


9.62 ± 1.01 


9.86 ±3.12 


RXJ1347.5-1145 


0.451 


13.0 


116 


1821 


2.17 


0.60 


5.00 


5.93 


3.77 


1446 ± 100 


14.03 ± 0.69 


12.05 ± 0.82 


13.83 ±2.82 


MACSJ1621.3+3810. . 


0.461 


34.7 


263 


2059 


2.29 


0.67 


4.96 


16.73 


1.00 


992 ±71 


7.53 ± 0.41 


6.42 ± 0.48 


4.37 ± 1.01 


MACSJ0329.6-0211 .. 


0.450 


26.6 


50 


466 


2.14 


0.49 


1.59 


0.00 


1.00 


915 ± 83 


5.24 ± 0.38 


5.23 ± 0.60 


3.64 ± 1.02 



^ The obtained best-fit parameters of the 3D gas density distribution. For definition see Eq.3 in lVildilinin et alll200.5tl) 



1 0^ I \ 1 — I I M I I I 




10 102 IQ3 

r, kpc 



Fig. 1.— Observed emissivity profile for MACSJ0159.8-0849. The solid 
line shows the best-fit profile and the dashed lines correspond to its 68% CL 
uncertainties 



detected in the Southern part of the Chandra image and is 
undetectable in the Northern part. In this case, we excluded 
the data within the lower (Southern) 50% of the field of view. 
We also excluded from the analysis the two outermost bins in 
the MACSJ1621. 3+3810 temperature profile because of po- 
tential contamination by the foreground emission. In addi- 
tion, there is a small-scale extended source at RA = 245.296°; 
Dec - 38.209°, also masked out. The only cluster in our sam- 
ple that shows some signs of a recent merger is RXJ1347.5- 
1 145 (Allen et al. 2002; Gitti & Schindle r 2004). In this case, 
we masked out the X-ray subclump at 20" Southeast of the 
cluster center. We also excluded the sector PA = 180° - 225°. 
This sector shows an excess of the X-ray surface brightness 
compared to all other directions, which is probably associated 
with the trail of the merging subcluster. 

The three-dimensional gas density and temperature profiles 
were derived from the azimuthally averaged X-ray bright- 
ness and projected temperature profiles centered on the clus- 
ter X-ray brightness peak. The projected temperatures were 
measured using the spectra accumulated in the concentric 
annuli with ro„,/r,„ - 1.5. These s pectra were fit to 
the single-temperature MEKAL model (iMewe et all 119851) 



with Galactic absorption. The absorber's hydrogen column 
density was fixed at valu es provided by HI radio surveys 
iPickev & Lockmanlll990l) . The metallicity was allowed to 
be free in almost all annuli except the outermost ones, where 
it was fixed at an average value from the two outermost bins 
where direct measurements were possible. 

The X-ray brightness profiles were extracted from the im- 
ages in the 0.7-2 keV energy band in concentric annuli with 
foiit/ri„ - 1.1. Using the observed projected temperature and 
metallicity profiles, the raw Chandra counts were converted 
to the profiles o f the projected emission measure, J n^npdl, 
as described in Vikhlinin et alJ ll2005b,) . The profiles from 
different pointings were added using the statistically optimal 
weighting. The resulting EM profiles were used to derive the 
gas density profile. 

3. PROFILE MODELING AND DETERMINATION OF THE 
CLUSTER MASS 

The 3-dimensional gas density profile was derived by fit- 
tin g the observed emission measure profile to a model used 
bv lvikhlinin etai] (|2005b ) to describe the 3-D gas density dis- 
tribution in the low-redshift Chandra sample (see their equa- 
tion 3). This function is a modification of the so-called /3- 
model (ICavahere & Fusco-Femianoll978h . designed to inde- 
pendently fit the gas density slopes in the cluster center, out- 
skirts, and in the intermediate region, and also to allow for 
an additional emission component in the very central region. 
The statistical quality of the Chandra data for high-z clusters 
(an example is shown in FigHJ is sufficient to fit the same set 
of parameters that Vikhlinin et al. (2005b) used for the low-z 
clusters. The best-fit parameters of the gas density profiles are 
listed in Tabled The projected emission measure correspond- 
ing to the best-fit model in MACSJ0159.8-0849 is shown by 
solid line in Fig^ 

For derivation of the 3-dimensio nal gas temperature pro - 
file we also followed the procedure in'Vikhlinin et al.'('2005l|. 
Specifically, a model with a great functional freedom was as- 
sumed for the 3D profile (see equation 6 in Vik hlinin et al] 
i2005b) . It was projected along the line of sight using the 
best-fit gas density profile. The projection uses the method 
that accurately predicts the single-temperature fit for a mix- 
ture of plasma spe ctra with diff'erent T (Vikhlinin 2005^ 
iMazzot ta et al. 2004). The parameters of the 3D model were 
obtained by minimizing x'^ computed using the projected 
model and observed temperature values at each radius. In 
some cases, the functional freedom of the 3D model was in- 
sufficient to accurately describe the data simultaneously in the 
very central region and at large radii. In such cases, we ex- 
cluded several innermost temperature measurements since our 
prime goal is the mass determination at large radii. The ob- 
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Fig. 2. — Observed temperature profiles of a sample of distant clusters observed with Chandra. The solid lines show the best-fit 3-D temperature profiles and 
the dotted lines coiTespond to their 68% CL uncertainties. The dashed lines show the projected best-fit 3-D temperature profiles (2-D). For MACSJ0329.6-021 1, 
we do not show the two outermost temperature bins because they can be contaminated by the foreground emission (see text); these bins were not used in the 3D 
modeling. 



served temperature profiles are shown along with the best fit 
models in Fig|2 

The total cluster mass profile was derived by direct appli- 
cation of the hydrostatic equilibrium equation (e.g., Sarazin 
1988) to the best fit gas density and temperature profiles. The 
mass profile was used to find the radius corresponding to the 
mean spherical overdensity A = 500 relative to the critical 
density at the cluster redshifts. The corresponding masses, 
M500, are reported in Tableland the values of rsoo are indi- 
cated in Fig^andl^] 

Using the best-fit three dimensional profiles, we also com- 
puted the average temperatures excluding the central region, 
as was done in Vikhlinin et al. (2005b) for the low-redshift 
clusters. We derived the average spectroscopic temperature, 
Tspec, a value obtained from a single-temperature fit to the in- 
tegrated spectrum, and computed the gas mass-weighted aver- 
age r,„g, both in the 70 kpc-rsoo radial range. These averages 
were used by Vikhlinin et al. ( 2005b) to derive the M-T rela- 
tions for nearby clusters and therefore we can use their results 
as the low-redshift reference. 

All the essential parameters derived from the 3-dimensional 
modeling are reported in Table[l] The parameter uncertainties 
were estimated by the Mo nte-Carlo simulations te chnique de- 
scribed in § 3.3 and 3.4 of lVikhlinin elal] J2005bl) . 

4. M-r RELATION 



TABLE 2 

Power Law Fit to Combined Chandra and XMM-Newton 
Mass-Temperature Relation 



Temperature average 


a 


r 


Tspec 


1.02 ±0.20 
1.33 ±0.20 


1.55 ±0.14 
1.65 ±0.15 



Note. — Parameters y characterize the slope of the M-T relation, £'(;;)M5oo ^ , 
and evolution in its normalization, Msqq/T^^^ oc E{z)~" . See text for details. 



Figure[5]shows the derived M- T relatio n for z > 0.4, which 
combines the XMM-Newton results from iKotov & VikhUninI 
(■2005) and the Chandra results described here. There is 
clearly a good overall agreement between the two datasets. 
There is a small difference in our modeling procedure for 
the XMM-Newton and Chandra data (we used more restric- 
tive models for Pg{r) and T(r) to fit the XMM-Newton data). 
However, we checked that this difference results in a negli- 
gibly small offset in the mass estimates. Applying the more 
restrictive models that were used in the XMM-Newton analy- 
sis to the Chandra data, we obtain that the derived values of 
Msoo change by 3% on average. 

To quantify the high-redshift M-T relation, we model 
separately the evolution in normalization and the slope of 
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T(keV) 



lower-right corner) at rsoo for the combined XMM-Newton(open circles) and 
Chandra(solid circles) sample. The solid lines show the low-redshift results 
from iVikhlinin et ah t2005h) . 



de-evolved relation. The evolution in the normalization 
is parametrized as M^m/T^^^ oc Eizy", where the low- 
redshift normalization is fixed at the Chandra measurements 
by VikJilinin et al. (2005b). The joint fit to the Chandra and 
XMM-Newton results for distant clusters gives a - 1 .02 + 0.20 
and a - 1.33 + 0.20 for the M - T^pec and M — T^g relations, 
respectively^ (see also Table|2}. These values are close to, and 
consistent with, the expected self-similar evolution, a - I. 

Assuming that the normalization of the M - T relation in- 
deed evolves self-similarly, we can "de-evolve" each mass 



measurement using the correction factor E{z) and fit the slope 
in the high-redshift M - T relation, E(z)M5oo oc T'*'. For 
the M - Tspec and M - T,„g relations we measure slopes of 
y - 1.55 + 0.14 and y - 1.65 + 0.15, respectively ( see also 
Table|2ll. These values are also close to the predictions of the 
self-si milar theory, y ^ 3/2, and to the low-redshi ft measure- 
ments dVikhUnin et alJ2005btlArnaud et alJ2005l) . 

5. CONCLUSIONS 

Deep Chandra and XMM-Newton observations of relaxed 
clusters at z > 0.4 provide direct hydrostatic mass measure- 
ments at the critical overdensity level A = 500. The resulting 
determination of the high-redshift M - T relation is compara- 
ble to the low-redshift studies both in terms of the sample size 
and mass measurement uncertainties for individual objects. 

We find that the M - T relation at z = 0.4 - 0.7 is con- 
sistent with the self-similar evolution of that for low-redshift 
clusters. Parametrizing the evolution of the normalization as 
we measure a close to 1 (Table. O. The 
de-evolved relation, E{z)M - T, has the power law slopes 
close to y = 1.5 (Table |2ll. This is an important confirma- 
tion of one of the key predictions of the cluster formation 
theory (Ekeetal. 1996; Brvan & Norman 1998) and also a 
welcome news for cosmological interpretation of the cluster 
number density measurements at z > 0. 

We should caution, however, that our results are based on 
the data for highly relaxed clusters, which represent only a 
fraction of the total population, especially at high redshifts. 
Extrapolation of the observed M - T relation on the entire 
population requires further studies with the help of numerical 
simulations and independent mass measurements from gravi- 
tational lensing. 
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Uncertainties include that for the normaliz ation of the low-z mass- 
temperature relations from iVikhlinin et ani2005hh . 
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